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Nonequilibrium wetting transitions are observed in Monte Carlo simulations of a kinetic spin system 
£S) , in the absence of a detailed balance condition with respect to an energy functional. A nonthermal 

model is proposed starting from a two-dimensional Ising spin lattice at zero temperature with two 
boundaries subject to opposing surface fields. Local spin excitations are only allowed by absorbing 
an energy quantum (photon) below a cutoff energy E c . Local spin relaxation takes place by emitting 
a photon which leaves the lattice. Using Monte Carlo simulation nonequilibrium critical wetting 
transitions are observed as well as nonequilibrium first-order wetting phenomena, respectively in 
the absence or presence of absorbing states of the spin system. The transitions are identified from 
the behavior of the probability distribution of a suitably chosen order parameter that was proven 
^ C*| useful for studying wetting in the (thermal) Ising model. 
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I. INTRODUCTION AND MOTIVATION 



Nonequilibrium wetting has been the subject of profound investigations in the past few decades [lH5| . In many of the 
models studied, the growth, or depinning, of an interface is described relative to a (usually one-dimensional) substrate. 
Temporal and spatial correlations in the interface are examined and dynamical universality classes are identified. In 
many cases, the nonequilibrium character of the phenomenon can be related to the breaking of detailed balance of 
configurational moves. In this paper we are concerned with one of the simplest ways in which detailed balance can 
be broken, leading to an intrinsically nonequilibrium system. In particular, after a move in configuration space the 
system can get trapped in certain configurations when the probability for the reversed move is identically zero. In 
thermal equilibrium, at finite temperature T, the reversed move always has a nonzero probability, proportional to 
O . the Boltzmann factor, which features the (finite) energy difference of the initial and final configurations. However, if 
we leave thermal equilibrium, by imposing constraints on the local absorption or emission of energy, detailed balance 
■ may be broken. We can go one step further along this line and leave the thermal context altogether by considering a 
^ ' classical system, say, at zero temperature, and providing a nonthermal mechanism for local energy exchange. 

To concretize our proposal, consider a (quasi-)two-dimensional lattice spin system at T — which is exposed to 
a photon bombardment from some external source. The photon energies hu are limited by a cutoff hv max = E c . 
We assume that a spin hit by a photon may absorb an energy E < hv, so that in all cases E < E c . Conversely, a 
spin may (always) relax by emitting a photon of arbitrary energy and we assume that photon leaves the plane so 
that the probability for absorption of emitted photons is negligible. The origin of the energy cutoff in this model is 
quantum mechanical. Although it is not necessary to invoke quantum mechanics explicitly to provide nonthermal 
energy fluctuations (random-field or random-bond disorder, electromagnetic fields, mechanical or chemical oscillators 
being alternative sources), it is a convenient frame- work for obtaining a sharp energy cutoff. In this manner we arrive 
at a model in which excitations of energy superior to E c are excluded, which implies that certain configurations can 
k>( | be trapping or "absorbing" . 

In the following we develop this model further and investigate how the character of a wetting transition is modified 
when thermal fluctuations are replaced by constrained nonthermal ones. We do so using Monte Carlo simulation and 
start within the context of the exactly solved wetting transition of a system in thermal equilibrium. Our paper is 
structured as follows. In Section II we test our simulation approach on the critical wetting transition in the two- 
dimensional Ising model @ . Section III is devoted to the definition of the nonthermal model, the analysis of the bulk 
phases and the observation and characterization of nonequilibrium wetting transitions of various nature. Conclusions 
arc drawn in Section IV. 



II. WETTING TRANSITION IN THE 2D THERMAL ISING MODEL 



Consider the two-dimensional square lattice Ising model with ferromagnetic nearest-neighbor interaction J > 0, at 
a temperature T below the bulk critical temperature T c and in zero bulk magnetic field. In this situation the bulk 
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consists of large coexisting regions of positive and negative magnetization. In the thermodynamic limit, the behavior 
of the bulk is independent of the boundary conditions, but this is not the case for the interface between the coexisting 
phases. In the case of a wetting transition, the surface excess free energy depends in a singular way on a surface field. 
To get efficient computational access to this transition, we use the same setup as in Q: a two-dimensional L\ x L-2 
square lattice of spins, periodic boundary conditions along the X-axis, open boundary conditions along the K-axis, and 
an anti- symmetric surface magnetic field Hi > acting on the spins along the open boundaries (see Figured]). This 
set-up is often referred to as one with "opposingboundaries" or "competing walls" and possesses surprisingly subtle 
and rich surface and bulk cooperative behavior [8l-ll2| . 

In the thermodynamic limit, letting L\ — ► oo followed by Li — ► oo, a sharp surface phase transition occurs as a 
function of the control parameter H\, assuming fixed T < T c . For small and opposing surface fields Hi and -Hi, with 
Hi below the wetting point H™ , the interface is localized at one of the boundaries and this two-fold degenerate state 
is called partial wetting: Figure[T]b). For Hi > H™ it is (free-)energetically favorable for the interface to wander away 
from the boundaries and the system is in the complete wetting state: Figure [1] a). For every temperature T <T C such 
a wetting point H™ (T) exists for which the excess surface free energy is singular and for two-dimensional equilibrium 
systems with short-range interactions this transition is known to be typically of second order. 




FIG. 1: Sketch of typical configurations during Monte Carlo simulations: dark spots are spin down, light ones are spin up. a) 
Large surface field Hi, complete wetting; b) Small surface field, partial wetting. 



Usually, to characterize accurately a phase transition in a simulation approach, one studies systems of different 
(large) sizes and analyzes the results adopting a finite-size scaling strategy. For our current reconnaissance study, 
however, it suffices to verify the order of the transition and to obtain an estimate of its location in parameter space. 
Therefore, a simpler approach is adopted. First, a bulk-like order parameter, the "magnetization amplitude" A is 
introduced 



A = — 



LlL2 



(1) 



where s(x,y) is the spin observable (= ±1) at site (x,y). For a system in the partial wetting state with an interface 
close to the boundary, A will be close to its maximum value of 1, while for complete wetting with an interface near 
the middle of the strip, A will be (much) smaller. Next, as a crude approximation to the successive limits L\ — >• 00 
and L2 — > 00 the fixed values Li = 100 and L2 = 10 are used, and simulations are executed for fixed temperature 
T and different values of H\. During each run a time averaged probability density P(A) of the order parameter is 
approximated by a normalized histogram. 

In Figure [5] the resulting histograms are shown for T/T c = 0.5. For Hi significantly smaller than 0.89 J the 
distributions are centered around a value close to unity, while for substantially larger surface fields, it is very unlikely 
to measure a A close to 1. Around Hi = 0.89J one observes a transition region where the distributions are broad, 
indicating the presence of large fluctuations consistently with a second-order interfacial phase transition. If we identify 
the wetting transition in the Ising model heuristically with the point where the distribution of A has a maximal 
variance, the result is in quite satisfactory agreement with the exact location. Note that in the limit T I the wetting 
transition, at Hi = J, is, exceptionally, of first order and purely determined by minimum energy considerations. In 
Figure[3]our estimates for the transition point for several temperatures are compared with the exact result of Abraham 
for the critical wetting phase boundary [6|, 



2J/feT = ( C0Sn ( 2 j/fcT) - cosh(2Hi/kT)) sinh(2 J/kT) (2) 



For temperatures not too close to T c we are able to determine the wetting transition fairly accurately. Closer to the 
critical temperature the interface is more fuzzy and it becomes increasingly difficult to differentiate between partial 
and complete wetting in a small system as we used. The method appears to be useful in principle for locating the 
wetting transition qualitatively. This is corroborated by the fact that the shape of the histogram provides an indication 
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FIG. 2: Histograms, or probability distributions, P(A) of the order parameter A, for the two-dimensional Ising model with 
Li = 100 and L2 = 10. Each curve is produced from a simulation of 5.10 s Monte Carlo steps/site. The suggested location of 
the wetting transition is at H\ = 0.89J, for which the variance of A is maximal. 




FIG. 3: Phase diagram of the wetting transition in the two-dimensional Ising model, in the variables reduced temperature 
T/T c and reduced surface field Hi/ J, where J is the bulk nearest-neighbor interaction. P: partial wetting phase, C: complete 
wetting phase. The full line is the exact solution |6(], the squares with error bars are the estimates from our Monte Carlo 
simulation and criterion of maximal order parameter variance. 



of the order of the transition. If, like in the thermal Ising case, one only observes distributions with one maximum 
and notices an increase of the variance in between the two sets of sharper shapes corresponding to the partial wetting 
(P) or complete wetting (C) states, one is in all likelihood dealing with a second-order phase transition. Close to a 
first-order transition, on the other hand, one rather expects distributions with two maxima that exchange dominance 
on crossing the transition point [l3j]. Further, the valley between the maxima is an expression of the hysteresis 
effect, the strength of which is proportional to the peak-to- valley ratio. Hysteresis is a clear signature of a first-order 
phenomenon. We will make use of these criteria in the analysis of the nonequilibrium model that we now introduce. 
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III. A NONEQUILIBRIUM WETTING TRANSITION 
A. Definition of the model 

In our preliminary simulation study of the (thermal) Ising model in the previous section the statistics was defined 
through a dynamics described by a continuous-time Markov process. In such an approach one aims to replace the 
canonical-ensemble average of the system by the time average in the long-time limit. Usually a dynamics is chosen 
that locally changes spin configurations {s} with rates that obey detailed balance with respect to the Gibbs measure 
defined by the energy functional E({s}). A physically appealing choice for realizing this is the single-spin flip Glauber 
dynamics [lij . One chooses at random a spin and calculates the energy difference AE of the system upon flipping 
the spin. Next, time is increased by one unit (usually 1 /volume) and the spin under consideration is flipped with 
probability Pm P - 




(3AE 



FIG. 4: Traditional Glauber spin-flip probability as a function of the reduced local energy difference of final and initial states, for 
the Ising model in thermal equilibrium. Note that the spin-flip probability is non-zero even for energetically highly unfavorable 
flips, consistently with the detailed balance property. 

where f3 = 1/kT, with k the Boltzmann constant. On repeating these spin-flip trials, starting from any initial 
configuration, in the long-time limit the system will visit all configurations with the correct Gibbs probability, i.e., 
the stationary state (probability distribution) of the Markov process equals the Gibbs measure. 

We now introduce a different spin-flip dynamics for which the stationary state cannot be described by a Gibbs 
measure for any energy functional. In such a system the dynamics is not induced by thermal effects and it is therefore 
an intrinsically nonequilibrium system. It is well known that nonequilibrium phase transitions can occur 

mm, ie -' 

expectation values of observables in the stationary state can depend in a singular way on control parameters. Our 
study focuses on a nonequilibrium transition of wetting type. 

From the probability function (|3]) it is clear that even if the energy cost for a spin flip is large, the flip probability 
remains strictly positive. This is an essential property of equilibrium models (not only of the Glauber dynamics) and 
it is related to the ergodicity of the dynamics. Therefore, a model without this property can be expected to have 
qualitatively different behavior. This leads us to the following definition of a nonequilibrium model. 

Consider the two-dimensional Ising model with nearest-neighbor interaction energy J on the square lattice. We 
use the same setup as in Figure [TJ but allow nearest-neighbor spins on the boundaries (y = 1 or y = L 2 ) to have a 
different nearest- neighbor interaction energy, Jr. The energy functional of this system can be written as 

E({s}) = -J ]T s(n)s(f 2 )- J T s (ri)s(r 2 ) 
(Pi,? 2 )gr <ri,f 2 )er 

-H 1 J2s{x,l)+H 1 J2s(x,L 2 ) (4) 

X X 

where J and Jr are the nearest-neighbor bulk and surface interaction energies, respectively, and H\ > is the (anti- 
symmetric) surface field. The two boundaries are denoted by T and located at y — 1 and y = L 2 . On this model we 
define the dynamics, which is not induced by thermal effects but by local absorption or emission of energy quanta. 
As in the Glauber model, we allow only a single spin flip at a given time. This flip can be the result of a photon-spin 
collision (absorption) or a spontaneous photon emission. We assume that the absorption probability decreases with 
energy, mimicking a certain photon frequency distribution, and becomes zero at a finite energy cutoff E c , which is the 
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maximum photon energy. We further assume that relaxation of a spin towards a lower local energy yields a photon 
that is emitted out of the lattice plane, so that we ignore spatial correlations between spin flips and avoid possible 
(secondary) absorption of energies greater than the cutoff. 

The single spin-flip probability function we propose is given by 



P aip = - 



1 + sign 



AE \ 
' E C J 



AE 



E c 



■1 



(5) 




AE / E n 



FIG. 5: Proposed spin-flip probability as a function of the reduced local energy difference of final and initial states, for the 
nonthermal Ising model. Note that the spin- flip probability is strictly zero when AE exceeds the energy cutoff E c , violating 
the detailed balance principle for thermal equilibrium. 



The important difference between our spin-flip function (Fig. 5) and the traditional one (Fig. 4) is not so much in 
the form of the curve, which we have chosen piecewise linear just for the sake of simplicity, but concerns rather 
the presence of a cutoff energy for excitations beyond which no spin flips can occur, since Pfa p (AE/E c ) = for 
AE > E c . This has major implications for the possible transitions between configurations. It can cause the system 
to be nonergodic and, as we will show, absorbing states can occur. Qualitatively, one can ask whether the energy E c 
plays a role similar to that of the energy 2kT in the thermal model, by examining the spin-flip functions and their 
slopes at AE = 0. On the one hand it seems irrelevant whether we call the energy scale, relative to which we measure 
the interactions, 2kT or E c . When both are large compared to the interaction and field energies the spin-flip function 
is sampled only over a small domain near the origin, where its value is about 1/2. In this regime our model cannot 
behave differently from the high-temperature limit of the Ising model. In the opposite limit, when kT and E c are 
both small compared to J, Jr and H±, however, the two models are physically different. In this limit excitations are 
excluded in the nonthermal model, while they are not in the low-temperature regime of the Ising model, although 
their probability is small. (The extreme case T — displays some special effects, which will be discussed separately 
further on.) While minimal energy considerations suffice, in the thermal model, to discuss phase transitions at low 
T, there is no such thing as an energy minimization principle in the nonthermal model for small E c , because the 
dynamics takes the system to one out of several absorbing states, regardless of their energy and only depending on 
which one is encountered first. 

Even greater physical differences between the two models arise when stable (i.e., "attracting") non-absorbing states 
and absorbing states are cohabitant (rather than "coexistent" , a term which has a well-defined meaning in the context 
of equilibrium phases) , as for example in the contact process at sufficiently large infection rate fl5l | . This regime is 
found at intermediate values of E c , for which ratios like J/E c are of order unity. Under these circumstances there is 
little or no connection between the behavior of the thermal and nonthermal models, in addition to the fact that free 
energy considerations do not apply to the latter. It is in this regime that our search will be conducted. 

Before we turn to a first analysis of the model, it is convenient to define dimensionless variables for the couplings 
and the reduced surface field: 



E' 



h ' = E:; 



(6) 



Analysis of the bulk phases 



Since we are interested in wetting transitions, in which an interface can form and move, the bulk of the system 
needs to be in the (two-fold degenerate) ordered phase. Therefore we start by analyzing the bulk behavior of an 
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infinitely large system, defined by the coupling parameter K = J/E c . It is instructive to recall that the equilibrium 
critical point of the thermal Ising model is located at J/kT c rs 0.4407 (square lattice). We have implicitly reproduced 
this in our simulations of the critical wetting phase boundary (Fig. 3), which terminates at T = T c . In the disordered 
phase limit, kT ^> J or E c ^> J, the correspondence E c w 2kT holds, as we discussed. In the spirit of a "low-if" 
approximation, that is, assuming that the physics found in the high cutoff energy limit can be extrapolated, one would 
thus expect a dynamical critical point at K c w 0.22 for the nonthermal model. However, the appearance of absorbing 
ordered states above a certain value of K may drastically alter this guess, especially if absorbing states occur already 
for K < 0.22. In this case, ordered bulk phases are being favored by the dynamics and we may expect the dynamical 
K c to decrease. 

To examine this, in the presence of the cutoff energy E c , it is useful to identify first the type of spin with the 
highest excitation energy cost AE for flipping: a spin aligned with all its neighbors. If this cost, AE = 8J, exceeds 
E c the spin flip is prohibited. Consequently, if K > 1/8 the system has two absorbing states in which the dynamics 
is frozen, {7} and {J.}: a configuration with all spins up or all spins down. For K < 1/8 there are no absorbing states 
and we find ourselves in a situation akin to that of the thermal Ising model. Therefore, noncquilibrium effects may 
be expected to drive the dynamical K c downwards towards the value 0.125, but not lower. 

The next important threshold for K is the value above which a spin can also not flip if it is aligned with all but one 
of its neighbors: K > 1/4. Above this coupling strength infinitely many (L2 — 1 for a finite system) configurations 
become absorbing; any configuration with one horizontal interface separating a spin up from a spin down region is 
frozen. This situation is reminiscent of that of the zero-temperature Ising model, for which the complete wetting 
state is {Li — l)-fold degenerate and the interface "does not move". Note that for K < 1/4 an interface between 
up and down domains is not frozen. Consequently, in order to allow interesting interface dynamics and the presence 
of absorbing states, we turn our attention to the regime 1/8 < K < 1/4. In this intermediate-coupling regime the 
behavior of our model is likely to be qualitatively different from both the finite-temperature and the zero-temperature 
Ising model. 

Now, for consistency, we need to verify that the dynamical critical value K c above which the system is in an ordered 
phase is low enough for the regime of bulk order to overlap substantially with the interval 1/8 < K < 1/4. Since for 
K > 1/8 the ordered states {7} and {1} are absorbing, it seems that K c > 1/8 is a plausible lower bound. In the 
absence of ergodicity it is not a priori clear that a random initial state at 1/8 < K will evolve into one of the two 
absorbing states. It is possible that a third, disordered, stationary state exists which is a dynamical attractor and 
features an order parameter probability distribution in which the probability to find the system in an absorbing state 
is zero. 

To decide on order versus disorder in bulk, we performed a first simulation of the nonthermal model, with uniform 
couplings K — Ky and h\ = 0. We used a square with equal sides of length L and imposed open boundary conditions 
to prevent the system from getting stuck in one of the states {7} or {1} (for open boundary conditions border spins 
have only three neighbors and can still flip against all of them as long as K < 1/6). As a simple and efficient order 
parameter we used the density of broken bonds pbr- 



with 5 the Kronecker delta, so that a broken bond is, as usual, defined as a pair of neighboring anti-aligned spins. 
In the infinite system limit, pbr is zero in the ordered phase and strictly positive in the disordered one. The scaled 
variance of this observable, i.e., multiplied by the number of spins L 2 , is akin to a dynamical version of the specific 
heat (capacity) per spin. If the bulk transition is of second order, we expect this dynamical specific heat to diverge 
approaching the bulk dynamical critical point. For finite systems we expect a peak, which becomes more pronounced 
and moves closer to the correct critical point with increasing system size. 

In Figure |5] the result of the simulation for three system sizes is shown. Besides the weak singularities in K = 1/8 = 
0.125 and K = 1/6 — 1.66..., which reflect effects of the sharp energy cutoff E c , one clearly observes the building 
up of a divergence, suggesting a second-order transition with 0.125 < K c < 0.13. To verify this and to get a precise 
estimate of the dynamical critical value K c , larger system sizes and an extrapolation are needed. Both ambitions are 
outside the scope of this paper. Nevertheless, from the previous considerations and these data we can convincingly 
infer that the X-window suitable for our explorations is: 



In the remainder of the paper we fix the bulk coupling therefore to a typical value in this range, K = 0.2, and we now 
embark on the investigation of possible wetting phenomena. Note that the spin-flip possibilities for spins on the two 
boundaries depend on the strengths of surface coupling and reduced surface field, and will be essential ingredients in 
our study. 




(7) 



0.13 < < 1/4. 



(8) 
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FIG. 6: Scaled variance of the density of broken bonds for the nonthermal model on a square lattice of linear size L, as a 
function of the bulk coupling parameter K. 



C. The wetting transition 

We return to the strip geometry (Fig.l). In order to understand the effect of the surface interaction energy Jr > 
and surface field Hi > it is, again, instructive to start with considerations involving the cutoff energy E c . The energy 
cost for a surface spin, aligned with all its three neighbors, to flip is AE = 4Jr + 2J — 2Hi or AE = 4Jr + 2J + 2H\, 
depending on the spin orientation with respect to the surface field. If both exceed E Cl the spins on both boundaries 
cannot destroy the absorbing character of the bulk and the two states {t} and {J.} remain absorbing. These states 
can be considered as extreme realizations of the partial wetting state, hence when 

4K r + 2K - 2hi > 1 (9) 

we concisely say that "partial wetting is absorbing" . For fixed K partial wetting becomes absorbing if the (reduced) 
surface field becomes small relative to the surface coupling. However, this does not imply that partial wetting is also 
an attractor of the dynamics. 

The converse can also take place. Consider a spin state in which all the spins on the boundaries are aligned and 
parallel to the local surface field. This includes, e.g., the state sketched in Figure QJi) corresponding to a complete 
wetting configuration. We simply refer to all these states as "complete wetting" , regardless of whether the interface 
is in the middle of the strip or close to a boundary. The important point is that the interface is separated from 
each boundary by at least one row of spins (cf. the multiply degenerate complete wetting state of the T = Ising 
model). Now, when the reduced surface field h\ and/or surface coupling Kr becomes large enough relative to the 
bulk coupling K, the boundary spins are unable to flip even when the bulk is in the opposite spin state and the 
interface can never quite touch the border, but it can still can wander between the boundaries (since K < 1/4). The 
corresponding energy evaluation gives that "complete wetting is absorbing" if 

AK T -2K + 2hi > 1. (10) 

For K = 0.2 the regions of absorbing complete and partial wetting are drawn in Figure[7] In the area of large surface 
coupling and/or (reduced) surface field, there is an ambivalent "phase" in which both partial and complete wetting 
are absorbing states. If the system is in either of these states it can never get out again. The dynamics in this phase is 
very different from that in an equilibrium system and it is not possible to identify partial or complete wetting phases. 
A certain fraction of initial configurations evolve into the former and the remainder into the latter. To attempt to 
identify "stable" phases based on minimal energy considerations is totally meaningless in this dynamical system [l7| . 
To construct a phase diagram for the wetting transition we performed simulations for K = 0.2 and analyzed them 
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FIG. 7: Phase diagram for the wetting transition in the nonequilibrium model for a fixed bulk coupling K — 0.2. To the left 
(for small Kr) we encounter a second-order wetting transition from the partial wetting phase (P) to the complete wetting 
phase (C). The thick solid line connecting the "measured" points is a guide to the eye. For larger Kr the transition changes 
into a first-order phenomenon, defined by a hysteresis region between two limiting lines (thick dashed lines). In the triangle 
(C&zPa) complete wetting is an attractor and partial wetting is absorbing, which leads to "cohabitation" of dynamical states. 
Another cohabitation is found for larger Kr, still in the hysteresis region, where both complete wetting and partial wetting are 
absorbing states (G'a&iPa)- The thin solid lines denote the exact location of the limits of occurrence of absorbing states. 



using the histogram method described above for the Ising model. For fixed Kr we made scans along h\. In Figure [8] 
the result of such a scan is shown for Kr = 0.18. Qualitatively, the behavior is similar to that of the thermal Ising 
model at critical wetting, i.e., the transition from partial to complete wetting upon increasing h\ appears to be of 
second order. In the inset of the plot, one can see that the variance of the order parameter A (defined in (JTJ) has a 
maximum around h\ — 0.14, which we identify with the transition point. The same type of transition was found for 
arbitrary values of Kr below approximately 0.22. 

For Kr above 0.22 the wetting transition enters the area of absorbing partial wetting and its character is drastically 
changed. There is now an entire subspace in parameter space characterized by the existence of two cohabitant states 
(see Fig. [7]). In this region complete wetting can be an attractor of the dynamics while partial wetting is absorbing. 
Let us examine this now in more detail. 

Suppose the system is initially in the complete wetting state in the region labeled C&P a in Fig. 7. We find that it 
remains in the complete wetting state for a (very) long time. Now h\ is slowly decreased. Below some value of h\ the 
complete wetting state rapidly evolves into the partial wetting state. This happens in the region labeled P a in Fig. 7. 
Since partial wetting is absorbing the system is trapped in one of the absorbing states {f} or {J,}. In this transition 
a histogram of the order parameter A abruptly changes from a distribution with a positive variance to a Dirac delta 
distribution on A = 1. 

The impression that the transition we just discussed is part of a dynamical first-order phenomenon is corroborated 
by the presence of a hysteresis effect: if we now reverse our path and increase hi so that we re-enter the region labeled 
C&P a in Fig. 7, the system stays in the absorbing partial wetting state. Upon further increase of hi, when crossing 
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FIG. 8: Histograms of the order parameter A, for the nonequilibrium model with L\ = 100 and L2 = 10. Each curve is 
produced from a simulation of 2.10 6 Monte Carlo steps/site. Inset: variance of A for the different values of the reduced surface 
field hi. 



the boundary 4Kr + 2K — 2h\ = 1 between the regions labeled Ck,P a and C in Fig. 7, it evolves abruptly into the 
complete wetting state. Also this reverse transition is part of the dynamical first-order phenomenon. 

When one studies, in a thermal system, a first-order transition in a simulation, one can use the histogram approach. 
Around the transition point the system can be in two stable states. For a finite system, this means it will stay for 
some time in one of these states but every now and then it will jump into the other. In this way, one long simulation 
can probe the whole configuration space and one can construct a histogram that in the neighborhood of the transition 
will show two peaks representing the two stable states. On setting the control parameters further away from the 
transition values, one of the peaks will start dominating until the second peak disappears: this sets the borders of the 
hysteresis region. 

However, our nonthermal system is equipped with absorbing states from which the spin configuration cannot escape, 
even for small system sizes. Under these circumstances the simulation cannot probe the configuration space properly 
and a naive use of the histogram method fails. To overcome this technical problem we introduced random spin flips: 
with a very low probability P r we flip spins, independently of the energy cost. We have to be careful with this, since 
it weakens the action of the cutoff energy and might restore the second-order character of the wetting transition. If, 
however, for very small P r we find non-second order behavior, we interpret this as evidence that the action of the 
cutoff is physically responsible for it and that a first-order dynamical phenomenon takes place. 

In our simulation we used P r = 0.001 and a system of small size, with L\ = 20 and Li — 5. (The transition 
time between two cohabitant states grows very rapidly with system size and with 1/P r .) The results for Kr = 0.24 
are presented in Figure |H1 For hi > 0.18, in the region labeled C in Fig. 7, we found a single-peaked histogram, 
representing a complete wetting state. From the moment one enters the absorbing partial wetting region labeled 
CSzP a , at hi — 0.18, a Dirac delta peak at A = 1 appears next to the peak of complete wetting. On decreasing 
hi this second maximum eats more and more of the probability distribution until at hi w 0.15, when entering the 
region labeled P a , the peak of complete wetting disappears. The values hi — 0.18 and hi w 0.15 are the two limits of 
the first-order transition phenomenon at the given value of Kr and mark the boundaries of the dynamical hysteresis 
region. These boundaries are indicated in Figure [7] by the dashed lines. We stress that the first-order phenomenon 
we observe pertains to an entire hysteretic region in the phase diagram, delimited by the two abrupt transitions we 
discussed. In our opinion, it is meaningless to attempt to locate a first-order transition "line" in the phase diagram, 
like one is used to do for equilibrium phase transitions. Based on our understanding of the dynamical behavior so far, 
we believe that a line of this sort does not and should not exist in our model. 

Figure [TU1 shows, for a system with parameters inside the part of the hysteresis region labeled C&P a , the time 
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FIG. 9: Histograms of the order parameter A, for the non-equilibrium model with L\ — 20, 1/2 = 5 and random spin flip 
probability P r = 0.001. Every curve is produced from a simulation of 2.10 6 Monte Carlo steps/site. 
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FIG. 10: The order parameter A as a function of time for the model with hi — 0.16, corresponding to one of the cases (middle 
right) for which -P(A) is shown in Fig. [5] 



evolution of the control parameter from which the histograms are constructed. One can clearly distinguish the two 
cohabitant states between which the systems jumps. Note that complete wetting is an attractor of the dynamics in 
the sense that even configurations close to the absorbing {f} or {1} states can evolve towards configurations with an 
interface in the middle of the strip. On the other hand, either one of the (quasi-)absorbing states is also frequently 
visited. Due to the non- vanishing random spin flip probability P r = 0.001, these visits are not permanent. 

Finally we remark that the separatrices in the phase diagram of Figure [7] between phases P and P a , or between 
phases C and C , are due to the singular behavior of the spin flip probabilities (|5|) about the value AE = E c and are 
irrelevant to the wetting phenomena. For instance, the transition from P to P a is a one-dimensional surface transition 
(which is not a sharp phase transition). 

In closing this section we also show, for completeness and for comparison, the phase diagram of the first-order 
wetting phase transition in the zero-temperature Ising model, in Figure 11. The equilibrium transition at H\ = J is 
shown, together with the spinodal lines at ±2J-p/J + 1 = Hi/ J. With respect to the Glauber dynamics all states 
with an interface parallel to the boundaries and the fully ordered states are absorbing states in the region P &C a in 
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between the spinodals. Beyond these spinodals, only either the fully ordered states {t} and {i} (P a ) or the interface 
states are absorbing (C a ). 
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FIG. 11: (Trivial) wetting phase diagram of the zero-temperature Ising model, in the plane of surface to bulk coupling ratio 
Jr/J and surface field to bulk interaction energy ratio Hi/ J. Partial (P; fully ordered, 2-fold degenerate) and complete wetting 
(C; straight interface parallel to the boundaries, (L2 — l)-fold degenerate) phases are separated by a first-order phase transition 
at Hi = J. The partial (complete) wetting state is metastable up (down) till the upper (lower) spinodal (dashed line). P a and 
C a denote absorbing phases with respect to the single spin flip Glauber dynamics, which in this zero-temperature limit only 
allows single spin flips that lower the energy. 



IV. CONCLUSIONS 

In this paper we studied the wetting transition from partial wetting to complete wetting in a nonthermal two- 
dimensional ferromagnetic Ising-like lattice model defined by a stochastic dynamics. The geometry is a periodic 
rectangular strip with opposing walls. We used a simple Monte Carlo simulation technique which could reproduce 
the order of the wetting transition and approximate the known phase diagram for the case of critical wetting in 
the equilibrium 2D Ising model with nearest-neighbor coupling J/kT and reduced surface fields ±H\/kT. For this 
equilibrium test case we employed a single spin-flip dynamics which obeys detailed balance with respect to the Gibbs 
measure for the given energy functional. Next we introduced an energy functional with, in addition to J and Hi, a 
tunable surface spin interaction energy Jr (between nearest neighbors) and introduced a non-equilibrium dynamics 
based on a piece- wise linear single spin- flip probability function which features an energy cutoff for excitations, E c , 
beyond which no spin flips can occur. Spin flips concur with photon absorption or emission, and kT is negligible 
compared to all other interaction and field energies involved. 

Analyzing the wetting transition using the Monte Carlo simulation technique, we found a rich phase diagram in the 
space {H\/E c , Jr/E c ), for fixed nearest-neighbor coupling strength in the intermediate regime 1/8 < K — J/E c < 1/4. 
In this range of K the interface between oppositely magnetized domains can still wander freely, but single spin 
excitations in a uniform domain are suppressed by the cutoff. Absorbing states can occur, and partial wetting as 
well as complete wetting can be absorbing. A second-order wetting transition is found for weak surface coupling, 
when absorbing states are absent. This transition changes to a dynamical first-order phenomenon, with hysteresis, for 
strong surface coupling. In a special region of the phase diagram, a dynamical attractor corresponding to complete 
wetting is cohabitant with an absorbing partial wetting state. The first-order phenomenon is qualitatively new and its 
properties are largely due to the presence of the energy cutoff for excitations. Only some of its characteristics can be 
understood as remnants of the trivial first-order wetting transition in the zero-temperature Ising model. The cartoon 
shown in Fig. 12 summarizes the most important features of our wetting phase diagram. 
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FIG. 12: Tutorial sketch of (part of) the wetting phase diagram of the nonthermal 2D Ising model, at fixed (intermediate) 
bulk coupling K = J/E c , as a function of surface coupling Kr = Jr/E c and surface field hi = Hi/E c . From the point of view 
of wetting phenomena there are three qualitatively distinct regions. A partial wetting phase is found for small hi , a complete 
wetting phase occurs at large hi and a hysteresis region appears for large Kr- The wetting transition line is of second order, 
except in the hysteresis region, which is bounded by two lines of abrupt dynamical transitions and which represents a dynamical 
first-order phenomenon. In the shaded region (in grey) the two extreme partial wetting states are absorbing. 
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